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ABSTRACT 


The rational use of optimality criteria is investigated for a class of 
structural synthesis problems where materials, configuration and applied 
load conditions are specified, and the minimum weight design is to be deter- 
mined , This study seeks to explore the potential of hybrid methods of 
structural optimization for dealing with relatively large design problems 
involving practical complexity. The reduced basis concept in design space 
is used to decrease the number of generalized design variables dealt with 
by the mathematical programming algorithm. Optimality criteria methods 
for obtaining design vectors associated with displacement, system buckling 
and natural frequency constraints are presented. A stress ratio method is 
used to generate a basis design vector representing the stress constraints . 
The finite element displacement method is used as the basic structural 
analysis tool. 

The optimality criteria are first derived for a general case and then 
modified for each type of behavior constraint. From these optimality 
criteria, recursive redesign relations are obtained for multiple constraints 
of the same behavioral type. In order to achieve high efficiency, design 
variable linking and temporary deletion of noncritical constraints are 
employed. The need for actual structural analyses is reduced by using 
first order Taylor series expansions to explicitly approximate the depend- 
ence of stresses and displacements on reciprocal design variables. Com- 
puter programs are written implementing some of the methods developed. 
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Results for several examples of truss systems subject to stress, dis- 
placement and minimum size constraints are presented. An assessment of 
these results indicates the effectiveness of the hybrid method developed. 
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CHAPTER I 


INTRODUCTION 

During the past fifteen years, considerable pro- 
gress has been made in the area of automated optimum 
design of structural systems. Optimum structural de- 
sign methods have been studied intensively and signi- 
ficant contributions have been made by many investiga- 
tors. These research studies can be usefully classified 
into three main categories as follows: 

(1) Application of mathematical programming algorithms; 

(2) Application of recursive redesign formulas based on 
fully stressed design methods and/or discretized 
optimality criteria concepts? 

(3) Mixed or hybrid optimization methods. 

In the first category, structural design problems are 
stated as inequality constrained minimization problems 
and solved numerically using mathematical programming 
techniques such as linear, nonlinear and dynamic program- 
ming. The most attractive feature of this approach is 
its generality in the sense that a broad class of struc- 
tural optimization problems can, in principle, be treated 
in a unified manner. This approach, pioneered by Schmit[l] , 
has actually enjoyed considerable success in a wide 
range of practical design problems [2,3,4 and 5]. How- 
ever, by about 1970 it became apparent that the applica- 
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tion of this method, in combination with finite element 
structural analysis, to large scale structural optimiza- 
tion problems required inordinately large numbers of 
analyses and long run times to solve problems of 
only modest proportions [6,7], 

This situation led some investigators to abandon 
the mathematical programming approach and direct renewed 
effort toward implementing recursive redesign procedures 
based on fully stressed design concepts and discretized 
optimality criteria. One of the first research efforts 
to focus attention on the discretized optimality criteria 
approach was reported by Venkaya, Khot and Reddy in Ref. 
[8]. This early effort has in recent years been followed 
by several notable studies [$,10,11 and 12] which pursue 
this same basic line of investigation. These studies 
have shown that the optimality criteria approach, includ- 
ing fully stressed design concepts, is well suited to 
achieving high efficiency in appropriate specialized sit- 
uations where only one type of behavior constraint dom- 
inates the optimum design. Indeed, during the past few 
years a widely held viewpoint has been that while mathe- 
matical programming methods are well suited to compo- 
nent design optimization, they are not practical for 
dealing with large structural systems. This assessment 
of the state of the art is well illustrated by the mixed 
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optimization method for automated design of fuselage struc- 
tures reported in Ref. [13]. In this work a fully stressed 
approach is used to obtain a gross overall distribution 
of material while the detailed design of rings and stif- 
fened panels is carried out using a mathematical program- 
ming technique. 

Nevertheless, other investigators, still attracted 
to the mathematical programming approach by its gen- 
erality, focused their efforts on a quest for the ef- 
ficiency improvements within the structural synthesis 
context [14]. With the same goal in mind, another 
type of mixed optimization technique was presented by 
Pickett . The reduced basis concept in design space 
and its initial exploration was reported in Refs. [15] 
and [16]. This work provides the foundation for a new 
group of mixed optimization techniques that will be 
called hybrid methods of structural optimization. Using 
the reduced basis concept it is often possible to dras- 
tically reduce the number of independent generalized 
design variables needed to obtain a good upper bound 
approximation of the optimum design. This is of cru- 
cial importance since the computational effort required 
to solve a mathematical programming problem grows rapidly 
as the number of design variables increases. It now 
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appears that the reduced basis concept in design space 
is one of the most promising approximation concepts in 

synthesis, although there are many open ques- 
tions that will require further study from both a theore- 
, tical and a practical point of view. 

The study reported here may be classified as a 
mixed optimization method. The objective of this inves- 
tigation is to explore the potential of hybrid methods 
of structural optimization using basis design vectors 
generated by optimality criteria and stress ratio methods. 
Primary attention is focused on the methods to be used 
in generating the basis design vectors. As previously 
mentioned, optimality criteria methods are well suited 
to the problem if only one type of behavior constraint 
is involved. It follows that the optimality criteria 
approach can ideally be used to generate basis design 
vectors if they are generated for each type of behavior 
constraint separately. Therefore, the optimality criteria 
approach is employed in this study as a tool for generating 
basis design vectors. Consequently, primary effort is 
concentrated on the derivation of discretized optimality 
criteria and the development of rational redesign proce- 
dures based on them. Although the basic statements of 
discretized optimality criteria are quoted from pub- 
lished descriptions [9-12] , this work attempts to estab- 
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lish a general viewpoint with respect to optimality cri- 
teria and their rational usage. 


5 



CHAPTER II 


PROBLEM FORMULATION 
AND 

SOLUTION METHOD 
2.1 Formulation of Problem 

Structural synthesis problems are considered for 
a structure which is modeled as an assemblage of a 
number of discrete finite elements. It is assumed that 
the types of elements used to model the structure are 
restricted to those which have the following special 
properties : 

(1) element stiffness is proportional to its mass; 

(2) there is one representative measure of intensity 
of stress state in each element. 

Elements of this type include axial force elements, 
shear panels, or constant strain triangles. The material 
of each element and the configuration of the structure 
are specified beforehand. All loads are applied 
directly or work equivalently at the nodes of the 
assemblage. The cross sectional area for an axial force 
element, and the thickness for a shear panel or a con- 
stant strain triangular membrane are taken as the design 
variables, and the minimum weight design of the structure 
is determined subject to both behavior and geometric 
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constraints. The former includes allowable tensile, com- 
pressive or combined stress, upper and lower bounds on 
joint displacement, lower bounds on general buckling 
loads and natural frequencies, and the latter includes 
minimum and maximum element sizes. Arbitrary design 
variable linking is also taken into account. 

Objective function 

The total weight of the structure is taken as the 
objective function to be minimized and it is given by 

W = ZI p; L; A i (2-D 

where NM is the fixed total number of elements represent- 
ing the structure, is the specified mass density of 
the element i, is the given length for an axial force 
element and the preassigned surface area for a shear 
panel or constant strain triangle element i. The A^ 
denotes the single design variable sizing each element 
i, namely the cross sectional area for axial force ele- 
ments or the thickness for shear panel and constant 
strain triangular elements. 

Constraints 

(1) Stress constraints 

For an axial force element the stress constraints 
are given by 
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where a . and a., denote the allowable compressive and 
ci ti 

tensile stresses for the element i, respectively, and 
ck j represents the stress in member i under load condi- 
tion j . For a planar element, shear panel or constant 
strain triangle, the stress constraint is given by 

O' .. < (2.3c 

~ C 

where o . . is the combined stress in element i under 
oi] 

load condition j, cr . denotes the yield stress of the 
material used in element i, and c is a safety factor 
(c > 1) . The combined stress is defined by 

i 

cr,, = c r/+ Ty - r* ay - jr/ y ) 2 (2.31 

where a and a are normal stresses in the x and y 
x y 

directions, respectively, and t represents the inplane 

xy 

shear stress referred to the x-y axes. 

(2) Displacement constraints 

Displacement constraints are defined by 


U L<J - i Uij i w ui j (2 - 4: 

where u T . . and u„. . are the lower and upper bounds on 
L 13 Ui] 

the ith joint displacement under the jth load condition, 
respectively, and u^^ is the ith joint displacement 
under load condition j. 

(3) Buckling constraints 

Any loading system applied on the structure can 
be expressed as 
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where Fj represents the load vector for the jth loading 

system, denotes the relative load vector in which 

the element of largest absolute magnitude has the value 

1 or -1, and oi. is a given scaling factor. Let p. 

3 *1} 3 

be a load vector corresponding to buckling of the struc- 
ture in the ith buckling mode under the jth load condi- 
tion. Then the buckling constraints associated with 
the jth load condition are given by 

2 c j (2.6) 

where c is a safety factor. At first glance, it may 
appear that the design procedure need only consider a 
single critical buckling mode. However, it has been 
shown that the buckling mode corresponding to the lowest 
load can shift from one mode to another during the 
design procedure. Therefore it is appropriate to 
include multiple buckling mode constraints for each 
load condition. 

(4) Natural frequency constraints 

Natural frequency constraints are given by 

(Jj- L £ (O) i (t)*u (2.7) 

where uk is the ith natural frequency of the structure, 
and w iL and w i(J are specified lower and upper bounds 
on the frequency, respectively. 
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(5) Size constraints 

Simple side constraints on the design variables are 
given by 

A; • < A- < A- (2-8) 

where A. . and A. are the minimum and maximum allow- 
lmin imax 

able element size for element i, respectively. 

(6) Linking of element sizes 

In actual design problems, any group of two or 
more elements may be required to: (1) have the same 

size, or (2) maintain fixed relative proportions through- 
out the design process. In order to implement this prac- 
tical feature, new design variables are introduced 
such that 

A; = OUj D; <2-9> 

Equation {2.9) means that element size A^ is controlled 

by the jth design variable D ^ , and its magnitude is 

given by a. .D.. It is to be understood that a- - is 

positive and that it is the linking coefficient between 

A . and D . . 
i J 

From equations (2.2) through (2.7), it is apparent 
that all the behavior constraints under consideration 
can be expressed in the following general form 

< ° ( 2 . 10 ) 

For example, a tensile stress constraint is given by 
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( 2 . 11 ) 


3* = <T; - <r ti £ 0 

where k is an arbitrary subscript used to identify a 
given constraint. Since the behavior constraints involve 
structural response quantities such as cr., u . p. 

l ij *i;j' 

and it follows that the g^, are functions of the 

independent design variables, that is 

?K = ?K ( 5 ) ( 2 . 12 ) 

where D represents the vector of independent design 
variables. 

Formulation of a design problem 

Finally the problem can be stated in the following 
general form: 


Minimize 


NM 

subject to 

w = 

S 1 
1=1 

h 4 

o , k 

= l. 2 

A i *1* 

S A i i 

A i m*/x 

A, - 

S' S ■ 



( 2 . 13 ) 


A.- - ft* Tv , { 

' 1 1 1 * I, 2. NPV 

where NC is the total number of behavior constraints, 
and NDV is the number of independent design variable 
Dj , which must be less than the number of elements NM. 
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2.2 Solution Method 

As shown in the preceeding section, the structural 
synthesis problem considered here can be expressed in 
the following mathematical form: 

Minimize W/ — W ( D ) 
subject to 

> S 0 , k = 1. 2, •••, NC (2.14) 

£ A i t..* , i ~ l. 2, ", 

9 m 

where D is a NDV dimensional design variable vector. 

The optimization problem represented by equation 
(2.14) can not be solved analytically, because for 
most practical structures the behavior constraints 
g^(D) cannot be expressed as explicit functions of the 
design variables. Therefore mathematical programming 
techniques are usually required to solve the problem. 
However, as mentioned in Chapter I, mathematical 
programming methods are not economically feasible for 
large scale problems in which the number of design 
variables exceeds one or two hundred. This is primarily 
due to the high dimensionality of the design space. 
Therefore, a simple and effective method to overcome 
this difficulty is to directly reduce the number of 
design variables. One of the promising methods for 
achieving this was presented by Pickett in Ref. [15], 
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and the essential idea is now employed here. 

Consider a set of basis design vectors D , £ = 1,2,... ,n 

Xj 

with n < NDV or n << NDV and define a set of generalized 
design variables 0^ such that 

D “ ZZ 0j D (2.15) 

Q«i * 

The optimization problem represented by equations (2.12) is 
now reduced to an n dimensional optimization problem in 
terms of the 0^, Jl = 1,2, ...,n. 

Minimize W = W C (f ) 
subject to 

iO, K I, 2, NC (2.16) 

* A^(0) — Av max , 1 - 2 , ***, MM 

As shown in Appendix B, the optimum design of the problem 
(2.16) is not, in general, the same as that of the problem 
(2.14) and it must be viewed as an approximation of the 
actual optimum design. 

To use the reduced dimensionality technique, a set 
of basis design vectors must be determined at the outset. 

They may be generated by a variety of methods. In this 
study effort is directed toward generating basis design 
vectors in the following way: 

(1) Classify the behavior constraints into several groups 
such as the stress group, the displacement group, 
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and so on. Then construct a set of subproblems , each 
of which has the same objective function as the 
original problem, while including only member size 
constraints and some of the behavior constraints. 

(2) Solve each subproblem by methods based on optimality 
criteria concepts as discussed in Chapter III, then 
take each subproblem solution as a basis design vector. 
After determining a set of basis vectors, there are 
several efficient mathematical programming methods that 
can be employed to solve the problem as stated in (2.16). 

In this investigation a program called CONMIN developed 
by Vanderplaats [17] , which is based on a modified feasi- 
ble directions method, was used. The overall optimiza- 
tion procedure employed herein is outlined in Fig. 1. 

The optimality criterion approach can be useful 
in its own right. However, these methods encounter 
practical difficulties when the optimum design is 
governed by more than one type of critical behavior 
constraint. Therefore, it should be recognized that 
optimality criteria methods may, in the future, be 
employed primarily as a source of basis design vectors 
for hybrid methods of structural optimization. 

On the other hand, as will be shown in the following 
chapter, problems involving size constraints and only 
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Fig. 1 Block Diagram of Solution Method 
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one type of behavior constraint can be treated effectively 
using the optimality criterion approach. In fact, includ- 
ing size constraints, in addition to a single type of 
behavior constraint, sometimes improves the computational 
efficiency. Furthermore, in those cases where a single 
type of behavior constraint dominates the optimum design 
of the complete problem, the appropriate basis design 
vector, generated including size constraints, will 
represent the actual optimum design exactly. Therefore, 
it is expected that the quality of basis design vectors 
may be improved if size constraints are included in 
each subproblem which generates a basis design vector. 
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CHAPTER III 


THEORY OF OPTIMALITY CRITERIA 
The purpose of this chapter is to formulate the 
optimality criteria principle and to develop the dis- 
cretized recursive procedure which will be used primarily 
to generate basis design vectors for each type of behavior 
constraint. The objective is to minimize the total weight 
of the structure, and all size constraints are simultan- 
eously considered in each case. 

In section 3.1 the basic optimality criteria concept 
and the associated recursive procedure is presented in a 
general manner. In the following secitons the same con- 
cepts are applied to each of several single behavior 
constraint types. 

3.1 General Concepts 

We will consider the following minimum weight design 
problem 

Minimize ***_*? 

W =• ZL fi Li A; (3.1a) 

\ 

subject to 

|k ^ ft > £ 0 , k - 1 , 2 , , NC (3.1b) 


A i min — A * — A *, m+.x ( 3 ,1c) 

A*' = a i j I>j , j = 1.2, , Nt>V (3. Id) 
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Using equation (3. Id), a design variable transforma- 
tion is carried out, leading to the problem statement 
expressed in terms of the independent design variables 
after linking, namely 

Minimize MDV 

W - H P; (3.2a) 

j-< 1 ‘ 

subject to 

J|< ( U> ) = 0 , |C = t. 2,-, Mt (3.2b) 


T> . ^ 

^ mt * — 


p i ~ 


(3.2c) 


where 


\ . = Mw f -A s mi*. 'N 


(3. 2d) 


TV,.... « Ml* ( Ad W0 v> "N 

V a... ) (3.2e) 

Now the optimality criteria for the problem given by 

equations (3.2a) through (3.2c) will be derived. Let 

D be a feasible design, i.e., 

J^S) 4 o , {or aflJL fc (3.3a) 


^ «ft.X , f 0>r } 


(3.3b) 


but not necessarily an optimal one. It is assumed that 
*?k is differentiable for all k and the following notation 
is adopted 


20 



(3.4) 



Consider a small change of the design, given by 6 d, then 
the corresponding change of the weight is given by 

HOV 

SW ~ §D; (3.5a) 

i * 

and the change of g^ is estimated by 


MDV 



It is possible to consider a number of small changes of 
design but among them we will consider only changes 
which do not violate any of the constraints. Such a 
change is called an admissible one and is defined as 
a change such that 

NDV 

$S>j i 0 , for k € k | ) 

SO, ?0, j •€ | j | - D,- J (3.6) 

SDj 4 0. f ^ ^ { j I t)j ® "D,'*.* i 

The design S can be improved if there exists an admissible 
change for which 6w is negative (see Fig. 2) . On the 
contrary, if the design is optimal, 5W must be nonnegative 
for all admissible changes. Figures 3 through 5 illustrate 
the typical cases of this situation schematically. Now 
the optimality criteria for the three cases shown in these 
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Fig. 2 Admissible Design Change reducing 
Objective Function 
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Fig. 4 Optimum Design, Case 2 
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Fig. 5 optimum Design, Case 3 
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figures are derived. 

First in Figure 3, the admissible direction in 
which the increase of the weight will be minimum is 
perpendicular to Vg^ that is 

3'-' SD . ^ 1,1 SD 2 = 0 (3 

Solving equation (3.7) for 6 D^ and substituting into 
equation (3.5a), yields 

SW = (-*/, + *0 SD* O 

N for ' 

In this case there is no restriction on 6 D 2 , therefore 
<5W >_ 0 for all possible 6 D 2 is satisfied if and only 
if the factor within the parentheses is equal to zero, 
i • © • / 


Aa/ | ^ — — • + 


0 


(3 


Let 




-A/7 


then from equations (3.9) and (3.10) 
Aa /~ 1 + ^ f I* > — 0 


(3 


(3 


A^z + * //,2 « 0 (3 

where X > 0 because from Fig. 3, w, > 0 and g. <0. 

X f JL 


. 7 ) 


. 8 ) 


9) 


10 ) 


11 a) 


lib) 
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In Fig. 4 , D 2 = D 2 min' therefore 6D 2 must be non- 
negative. In this case <5W > 0 is satisfied when the 
value in the parenthesis is not less than zero, i.e.. 


~AaTi 


h 


+ 4</~z Z 0 


(3.12) 


Using the same definition of A, see equation (3.10), it 
follows that 

4Af* + >?;,* i 0 <3 ' 13) 

From Fig. 5 it can be shown in a similar manner that 
5W > 0 is satisfied if 

+ A 4 0 (3.14) 

From the above discussion, the optimality criteria for 
these cases can be summarized as 

j - 0 / f* r \ * i I ^ <*>• 

+ * fl,j % - °/ V i 4 J’min (3.15a) 

v £ 0 , •j-**' - j t J" m*x 

and 

K > 0 (3.15b) 


For a general case, the optimality criteria can be derived 
on the same basis, that is. 


%VJ = 1 Z 4«/ySt>: S 0 

j" * 1 


for all 6Dj such that 


(3.16a) 
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MOV 

= ^ Jloj SD| 4 0 , for k Kact 




w i € 

J* win 

4 

0 , 

for j « 

vF 


(3.16b) 


According to Farkas ' lemma, a set of equations (3.16a,b) 
can be shown to be mathematically equivalent to the 
following (see Appendix A) , 


HC 

i; 


for 

j * J 

+ fa ^ 


0 , 

f° r 

j 4 3”mtn 

where 

< 

V * 

0, 

V 

fl ^ J" m*x 


(3.17) 



? 0 , for k <: K act 
“ 0 , for k K act 


Any design which does not satisfy (3.17) is not an 
optimum design because from (3.16a,b), there exists 
at least one admissible change of the design for which 
<5W < 0. In this sense the set of conditions stipulated by 
equation (3.17) is called the optimality criterion. 

A redesign equation can be obtained from the opti- 
mality criterion in various alternative ways [9], [10] 
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and [12], but all these methods are based on the follow- 
ing central idea. The optimality criterion given by 
(3.17) can be rewritten as 




’ = l , 


t € J 

N C 

J 

< \ 




- x / 



^ i. 


J sJ 


(3.18) 


Equation (3.18) suggests that if 1^ < 1 and D ^ > Dj m ^ n for 

a current design, must decrease to obtain an improved 

design at the next iteration, and conversely if 1^ >1 

and D. < D. , then D. must increase. On this basis, a 
3 D max ] - 

redesign rule is expressed as a function of I ^ . In 

this research effort the following plausible redesign rule 

is used 


V, 


-1 


(i, )* c r>j ) s , f°>- * o 

0 , 


(3.19) 


where ( ) represents the s-th cycle in the iterative 

s 

design process. In order to include the size constraints, 
the following relation must be appended to (3.19) 




Sf \ 


( D| )*+\ , i \ 

^ ( Dj- 

)st\ ^ £>j' m** 



( X>j ) ji. i ^ 

O j min 

(3.20) 
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The reason for using the square root of 1^ in equation 
(3.19) will be discussed in subsequent sections dealing 
with specific constraint types. From the foregoing dis- 
cussion it is obvious that if (D.) = (D.) for all j, 

3 s+1 3 s 

design (D) g satisfies the optimality criterion, and conse- 
quently it may be an optimal design. If (Dj) g+1 ^ (Dj) g 
for some j, the design can be improved further, and re- 
peated application of the redesign equation will converge 
toward a design satisfying the optimality criterion. 

The value of multipliers X^ can also be obtained 
from the optimality criterion. From the first equation 
in (3.18) 

mc a , 

1+ +•>- | « J (3.21) 

i 

Eliminating the terms corresponding to inactive constraints 
from equation (3.21), since those X k = 0, yields 

1 + 2_ =0, {or j J (3.22) 

- 

where k act represents k«-K act . Note that the values of 
Wj are known constants and assume that the values of 
the g^ j are known for all combinations of design var- 
iables j-ej and active constraints k€iK , then equation 
(3.22) represents a set of simultaneous linear equations 
in which the are the unknowns. Since the number of 
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design variables NDV is not necessarily equal to the 
number of active constraints, it is in general not 
possible to directly solve (3.22) for the A^» Conse- 
quently an indirect method must be employed to solve 

for the A, values, 
k 

First consider the following optimization problem 
in which the A^ are the unknowns 

Minimize X = 2 j i + Ate j (3.23) 

j -6 T fe »vt J 

If a solution is obtained which minimizes I and l m ^ n 
has zero value, then this solution gives a set of the 
optimal values for the multipliers A^, and the design 
cannot be further improved. If the minimum value of 
I is greater than zero, it means that the value in 
the parentheses is not equal to zero for at least one j, 
and therefore the design can be improved. The problem 
given by (3.23) can be transformed to the following 
linear programming problem 

NAt 

Minimize 2_ + J, ) 

Q. X= \ 


subject to 


H >*{ T. 


J 







20, XI 2 0 ^ ^ o 


(3.24) 
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Now this problem can be solved by the Simplex Method, 

a 

for example. The method to be used in estimating the 
value of the g k ^ will be discussed in subsequent sec- 
tions dealing with each of the particular constraint 
types considered. 

In executing the computation, however, special 
attention must be given to identifying the maximum or 
minimum size elements as well as active and inactive con- 
straints. First consider the problem of determining 
those elements that are to take on their maximum or min- 
imum values. Generally no information is available on this 
at the outset and therefore an iterative procedure is re- 
quired. The procedure is outlined as follows: 

(1) Initially assume that all design variable side 

constraints are inactive, i.e.,D. . <D. <D. 

jmm j jmax 

for all j, and compute the value of X^ using (3.24); 

(2) Calculate (D^) s+1 from equations (3.19) and (3.20), 
and use the results to identify the active design 
variable side constraints; 

(3) If the distinction remains unchanged, it is done, 
otherwise, use the new distinction to repeat the 
procedure . 

This iterative process must be carried out during each 
redesign step until the set of active side constraints 
has definitely stabilized and remains unchanged. 
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Next consider the problem of identifying the set 

of active behavior constraints . t The procedure discussed 

! 

in the foregoing paragraph will be applicable, however, 
it may waste a great deal of computational effort, because 
in structural optimization problems, a large number of 
behavior constraints usually need to be considered but 
a relatively small number of these constraints are 
active . 

The procedure is as follows: 

(1) Select several active constraint candidates if they 
are known; 

{2) Analyze the current design and evaluate all con- 
straints; 

(3) Find the most critical constraint, which may or 
may not be violated, and compare it with the 
list of preselected active constraint candidates; 

(4) If it is already an active candidate, continue 
with the current list, otherwise, add it to the 
list of active constraint candidates; 

(5) If the number of active constraint candidates exceeds 

the number of design variables, eliminate the one 

which is least critical during the redesign step. 

tThis is not necessarily required because even if inactive 
constraints participate in the process, they will be auto- 
matically eliminated by the computational result of ^ = 0 
for the corresponding k. 
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In the foregoing method, there is no need to compute the 
corresponding to inactive constraints, and this sub- 
stantially reduces the computational effort required. 

Now an entire optimization procedure based on 

optimality criteria concept is available and it is 

summarized as follows (see Fig. 6) ; 

(1) Select several active constraint candidates and 
pick an initial design; 

(2) Analyze the current design and evaluate all con- 
straints ; 

(3) Find the most critical constraint and determine the 
active constraint candidate group for the upcoming 
redesign step; 

(4) Compute the gradients of the active constraint candi- 
dates ; 

(5) Compute the value of multipliers A k from (3.24) and 
use the results to generate a new design from equa- 
tions (3.19) and (3.20); 

(6) Determine which elements are to take on their max- 
imum or minimum values and repeat (5) and (6) until 
the set of active design variable side constraints 
has stabilized; 

(7) Check to see if the new design satisfies the pre- 
scribed termination conditions. If so, go to (8), 
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Select N constraints as the 
active group (N NDV) . 


Pick initial design, 



Find the most critical 
constraint. 


'is the constraint 


included in the 
\active group? 



Analyze current design and 
evaluate all constraints. 


Add it to the 
active group and 
eliminate the 
least critical 
constraint from 
the active group. 



Add it to the active 
group . 


Compute the gradients of 
all active constraints. 


Generate new design. 


Classify the design 
variables as active 
or passive. 


Does the \ 

classification /no" 
change? / 


I — 1+1 


-> * 
I >1 



Compute the values of 


Does the design 
satisfy the 
termination 
criteria? 



Fig. 6 Block Diagram of Optimality Criteria Approach 


















otherwise, go to (2) and repeat; 

(8) Evaluate the remaining behavior constraints using 

the final design and stop. 

The purpose of step (8) is to check on whether or not 
the final design is the optimal solution of the whole 
problem under consideration and to obtain information 
for selecting active constraint candidates in the other 
subproblems. 

At this point some points which should be noted 
when applying the present method to an actual design 
problem are discussed. As is obvious from the deriva- 
tion, the optimality criteria are nothing more than 
the necessary conditions for local optimality, and they 
do not guarantee that the design obtained by the present 
method is the global optimum. An example of this situa- 
tion is shown in Fig. 7, where either D-^ D 2 or D* will 
be obtained. Among them, however, only D* is the opti- 
mum design, and the other two are obviously not optimal. 
But this situation can be avoided by setting up some 
limitations on the values of multipliers. At design 
D*, both constraints are active, and it follows that 
the values of A^ and must be non-negative. At design 
^ 1 • 9^ is active and g^ is inactive but violated, and at 
®2' ^1 i nact ive but violated. This situation suggests 
that the value of a multiplier A^ corresponding to a 
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Fig. 7 Illustration of Optimum Design and 
Infeasible Designs 
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violated constraint should be positive. Based on 
this fact, we set up the limitation on the value of 
multiplier as 

A* 2. e k ( > o ) (3.25) 

for violated constraints instead of ^ > 0 in (3.24). 

The value of £ k may be determined in various ways, but 
it will be better to define it as a function of (9^/9^*)'+ 
and the function should be defined for each problem. In 
the effort reported here 

88 ( -|jr - » ^ k 4 1 ?*> 1* \ < 3 - 26 > 

was used, where is an appropriate constant. 

Another example is shown in Fig. 8, where the opti- 
mality criteria are satisfied by design , E> 2 and § 3 . 
Among them is the global minimum, S 2 is a local mini- 
mum, and D 3 is a local maximum. An effective method 
of coping with this situation has not been found, and 
it would appear that this difficulty represents one 
of the current shortcomings of the optimality criteria 
approach. It should be noted, however, that the formida- 
ble difficulties posed by relative minima are not unique 
to the optimality criteria approach. 

tFor the notational convenience, we rewrite the constraint 
in the form g, - g, * <_ 0, where g * denotes the specified 
upper bound on g^. K 
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Illustration of Global Minimum, 
Local Minimum and Local Maximum 


Fig. 8 


3.2 Stress Constraint 


Even for the simplest class of structures , such as 
trusses, the optimality criteria approach developed in 
section 3.1 is not practical for stress constraints. 
However, there is a useful method, which is commonly 
used, that does not require the derivatives of constraints 
although it deviates somewhat from the general optimiza- 
tion theory. 

In structural design, it has often been assumed 
intuitively that the best design is one for which every 
mode of failure considered would occur simultaneously. 

From this idea it followed that for stress limited design 
problems the best design would be one in which each member 
is fully stressed under at least one load condition. How- 
ever » it was shown by Schmit [1] that the fully stressed 
design is not necessarily the minimum weight design. 
Nevertheless, the fully stressed design scheme still has 
practical significance because of the following charac- 
teristics: 

(1) The fully stressed design always coincides with 
the minimum weight design for statically deter- 
minate structures. 

(2) For the case of statically indeterminate structures, 
the fully stressed design may be a good approxima- 
tion that is often acceptable for practical purposes. 


40 



(3) The fully stressed design procedure is familiar 

and easy to apply in comparison with methods based 
on mathematical programming concepts. 

For these reasons the fully stressed design concept is 
adopted here as the source of design basis vectors for 
stress constraints. That is, "the optimum design for 
stress constraints is assumed to be one in which each 
member is fully stressed under at least one of the load 
conditions." 

The fully stressed design is usually obtained by 
the stress ratio method which is derived based on the 
assumption that the internal force distribution remains 
unchanged during modification of the design variables 
in each redesign step. This is equivalent to 

r Oa a! , ecu, 


C<Ti ), ( Ai ). 


t (T; c A* , ^ cot), < o 


(3.27) 


where A^* denotes the optimal design for . Conse- 
quently, the following redesign equation can be obtained 

<• Ai « C. (Ai), (3.28a) 


where 


C; 


IQ! ) s 


if cot- ) s S 6 


0 s ; c 


: f ), < o 


(3.28b) 
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(3.29a) 


For linked design variable D ^ , the redesign equation 
must be modified as follows: 

c x>j ) = Cj (Dj;. 

where 

C: = C; (3.29b) 

> U j 

Including size constraints, equation (3.20) must be con- 
sidered together with equations (3.29). The redesign 
procedure can now be summarized as follows: 

(1) Pick an initial design; 

(2) Analyze the current design and compute the stress 
ratio for each combination of elements and load 
conditions ; 

(3) Find the maximum stress ratio for each design variable; 

(4) Generate a new design using (3.29) and (3.20); 

(5) Check to see if the new design satisfies the pre- 
scribed termination conditions. If so, go to (6), 
otherwise, go to (2) and repeat; 

(6) Evaluate the remaining behavior constraints using 
the final design and stop. 

This procedure is shown in block diagram form in Figure 

9. 
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Fig. 9 Block Diagram of Stress Ratio Method 
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3.3 Displacement Constraint 

As shown in Chapter II, a displacement constraint 
is given by 

u Lij 4 u ;i 4 W uCj (3.30) 

In the majority of practical problems, is positive 

and u . . is negative. Assuming this is so, the con- 
straint can be expressed in the following form: 

* - . 


f T U " U, 


(3.31) 


where u is a displacement vector including u^ in its 
ith row, and ? is a unit force vector which has only 
one nonzero element in its ith row, namely 

t ' if U*'j 2 0 
. if U*j < 0 

since ? T u represents the absolute value of displacement 
and denotes 


f; = 


(3.32) 


* 


u <J *-j " ’ 


z 0 

U L ;. , i 

f 

< 0 


(3.33) 


Hereafter equation (3.31) will be taken as the displace- 
ment constraint form because this form facilitates the 
derivation of partial derivatives Og/9Dj) assuming the 
use of a displacement type finite element method of struc- 
tural analysis. For the sake of simplicity, a case with 
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only one constraint will be considered, let 

. "r*T -> * , 

then 


-24 - 

ao; >n- ^ ^ / 




J 

From equation (3. Id) 




- aii 




<1 pAi 


Since f is independent of the design variables 


_ ( V - \ ^ 

w > ’ f IaT 


M 


(3.34) 


(3.35) 


(3.36) 


(3.37) 


The static equilibrium equations for the displacement method 
of structural analysis may be written in matrix form as 

LK] u = F <3.38) 

where [K] is the stiffness matrix of the structure, and 
f is the external load vector. Assuming F is independent 
of and differentiating both sides of equation (3.38) 
gives 


3U 

9A.- 



aOO 

>A * 



(3.39) 


Substituting equation (3.39) into equation (3.37) yields 


rcr*)« -rw 




a A; 


u 


(3.40) 
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Now define a new vector such that 


r = [Kf ? 

then r represents the 
unit force vector ?. 

r T ■= ? t LK]" 


(3.41a) 

response of the structure to the 
Since [K] is symmetric 

(3.41b) 


For the class of structures considered in this study, 
the stiffness matrix can be expressed in the following 
form 


NH 

L K 1 « 2__ [ ki } A; (3.42) 

1 

where is the unit element stiffness matrix for element 

i, which is independent of the element size. Therefore 


a[K] 


[ k; 3 


(3.43) 


Substituting equations (3.41b) and (3.43) into equation 
(3.40) gives 

)= (3.44) 

From equations (3.35), (3.36) and (3.44), it follows that 



&*-i >* T L je; 3 U 


Using equation (3.17), the optimality criteria for 
single displacement constraint is obtained, namely 


(3.45) 


a 



Multiplying both sides of equation (3.46) by yields 
the following standard form 


W\ 


>U: * 0 , for -j-feJ 


where 


w i “ *'1 


(3.47a) 


(3.47b) 


u j = ?■ ft< J “ 


(3.47c) 


Now Wj represents the total weight of the elements in 
the group j, and IK is the internal virtual work in those 
elements associated with the jth design variable (D^). 
Equation (3.47a) can be rewritten as 


I 

X 




W 


( 3 . 47d) 


1 


Therefore the optimality criterion can be stated as "the 
ratio of the internal virtual work over the weight is invar- 
iant for all active element groups." 

For the case of multiple constraints, the criteria 

can be generalized as follows 

N< 


Wj - Xk Ui.j - 0, i < J 


k = I 


(3.48a) 


where 




= 0 , for k -S' K act 

= o , for k K a*-t 


(3.48b) 
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and U.. represents U. for the kth constraint. 

K J J 

Consideration is now given to the redesign equation. 
Initially the same assumption used in fully stressed design 
(see Section 3.2) is made, that is, the internal force 
distribution remains unchanged during modification of 
design variables in each redesign step. Then the follow- 
ing relations must be satisfied for the case of one con- 
straint, if 

(X Vs., = (3.49a) 

then 

\A/j ) S -M = C,' t Wj )s (3.49b) 

and 

( Uj) w - - 7 - (Uj), (3.49c) 

At iteration s, if W. - ALL ? 0 for some j which are 

assumed to be active, then (D.y . must be determined so 

3 s +1 

that 


- x(Uj> w - 0 

Substituting equations (3.49b), (3.49c) into equation 


(3.50) 


(3.50), and solving for , we get 




(3.51) 


If the value in the parentheses in equation (3.51) is 
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0. 


negative, it means that Dj must be inactive, i.e., Dj = 

For the case of multiple constraints 

\ 

Cj = Ij 2 (3.52a) 

where 

L i = 

Based on the foregoing discussion, the redesign equation 
can be summarized as 




= < 


Cj ( Dj 3 s , 


if Ij 5 o 
if Ij < 0 


(3.53) 


where Cj and Ij are defined respectively by equations 
(3.52a) and (3.52b). Including size constraints, equation 
(3.20) must be employed concurrently. 

Finally, consider the method used to compute the 
values of multipliers . Since the general idea was pre- 
sented in section 3.1, it is only necessary to discuss 
/ 

the method used to compute U, ./W . , which corresponds to 

the term of g, ./w. in section 3.1. For the active 

D 3 

constraints, the following relation must be satisfied: 

NDV 

LU; = ut (3.54a) 

i* 1 * 

NI>V 


( Ukj) = tUk), (3.54b) 

i - 1 1 
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where u k * is the specified value of the kth constraint 
and (u k ) denotes the corresponding displacement for 
the current design. They can be divided into active 
and passive parts as follows 

NOV 

u Nj = tL u itj + H u kj 


(3.55a) 


1 






H PV 


C < Utej ) s 2^ < U fej (3.55b) 

Again using the same assumption employed in fully stressed 
design 


(3.56) 


U h : = £ ( Uk) 

j -* 1 j-*T J £ 

Let u k ° denote the value of and let 

Uk i ’ ( U|, j ) *Tfe • *" r 

then satisfies equation (3.54a). Based on this fact 

the following equation can be used to estimate the value 

of U. ./W. , 

3 


(3.57) 


(m 


(3.58) 


'6 ^ u*p s < o 


because (U k j) g < 0 means that D ^ is inactive at least for 

the constraint. After obtaining the values of the U, ./W. 

*3 3 
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the A k are obtained by solving the linear programming 
problem defined by (3.24). 
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3.4 Buckling Constraint 

As indicated in Chapter II, a buckling constraint 
is given in the following form 

h = “ t* + f* £ 0 (3.59) 

and p^ is defined by 

iKltffc = f K [Ke]i7fc (3.60) 

where the subscript k denotes that the buckling load under 
consideration is the kth constraint. Let p., * denote 
the prescribed lower bound on p^, let u^ represent the 
corresponding buckling mode shape, and let [K„] denote 
the geometric stiffness matrix of the structure, which 
is symmetric and independent of element sizes for the 
class of structures considered here. 

Differentiating both sides of equation (3.60) with 
respect to A^, we get 

He d Uk - (3.6i) 

Premultiplying equation (3.60) by [auJ,/3A;] T and 
premultiplying equation (3.61) by Uj^, then subtracting 
the former from the latter yields 

[fci] U* = tf/ [K*] tT* (3.62) 

From equations (3.36), (3.59) and (3.62), we get 
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(3.63) 


= 


9Dj 


_ r Mtlkll 

;«j ^ < [K*] u* 


Substituting equation (3.63) into equation (3.17) yields 
Multiplying equation (3.64) by Dj , leads to the following 


optimality criterion, namely 
mc . 

At 


w r h 1u u >i = 0 ' ^ 


(3.65a) 


where 


Ufcj - 2_. u7 

^ i 

M* - u? [Ks] u„ 


(3.65b) 


(3.65c) 


The redesign equation can be obtained in exactly the same 
manner as for displacement constraints, and that is 

CjCDp, , ij I; so 


< D , - 


( 3 . 66 a) 


where 


, i-j- Ij < t) 


i- = ( r Ja.. 

<) \ fa Ms Wj y. 


(3.66b) 


and 


c i = V 


(3.66c) 
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Including size constraints, equation (3.20) must be 
employed concurrently with equation (3.66a). 

The value of Xy. can be estimated in the same manner 
too. For an active constraint 


NDV 


L- 




(3.67a) 


T (rtf). - < * >: 


(3.67b) 


They can be divided into active and passive parts as 
follows : 


* U tt j 2_» U fej + 2^ U Kj 


r 

£(u^.z: T (u k j). 

r' J a* J 


Again using the same assumption employed in fully stressed 


design 


u*j = 2_ (Ui.:), 

J j+T J 

Let p v ° denote the value of 7Z. (Uni )« , and let 
k )<7 } 


(3.69) 


. (M*L\ K- 

V Ml e A ( fn),- 


t A.® 


- 

r i* 0 


w j *■ j 


(3.70) 


then equation (3.67a) is satisfied. Based on this, it 
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follows that the following equation may be used to estimate 


the value of 


U>j 

MbW; 



A f* - K 
M«w; 4 t; 


if 


I»i SO 


_UKi_'v 
MnVv/j 's 




(3.71a) 


where 



c 


j . A 

M* W j Js 


(3.71b) 


55 



3.5 Natural Frequency Constraint 


A lower limit natural frequency constraint is 
given by 

J* = - i 0 (3.72) 

and q k is defined by 

[K] U> = % k [M] Ujr U.73) 

where the subscript k denotes that the frequency under 

consideration is the kth constraint. Let q, * denote 

^k 

the prescribed lower bound on q^, and let u^ represent 
the corresponding natural mode shape. The mass matrix 
of the structure considered is represented by [M] and 

for the class of structures considered here 

NM 

[M 3 - Zl A; (3-74) 

;=i 

where the fm^] are unit element mass matrices independent 
of the element size. 

Through a development that runs parallel to that 
used in the case of a buckling constraint, the following 
optimality criterion can be obtained 
c ^ 

Wj “ 2_- ~~ ( Unj “ SfcTkj ) « 0, f J (3.75a) 

where 




(3.75b) 
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(3.75c) 


T * = uj [m] U* ( 3 . 75d) 

The redesign equation is also obtained in a manner anal- 
ogous to that previously employed in the case of a buckling 
constraint and the result is 




Cj CDj) s , if So 
« , ij Ij < o 


where 


fe-i \ T* W'j 


(3.76a) 


(3.76b) 


and 

Cj - I- 1 (3.76c) 

If size constraints are imposed, then equation (3.20) 
must be used in conjunction with equations (3.76). 

For an active constraint, the following relation 
must be satisfied 


NDV 

2- ( Uk.- - ll TkI ) - 0 (3.77) 

i" 

Consequently 


( Ukj Ik] ) 2_ ( U i*j ) (3.78) 
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Now we assume that 


U K j “ CUhj), 'j 

T»j = ( T,J ) . i 




then 


.11 c u„ r t: Ti,j ) = - ^ c u»j - ?; T«j > s 


' 1+7 

Define the following two quantities 


Q 


~ ^ T| *i 


Q. p — 


- 2_ ( U»i - ti T '>j ), 

j*T J J 


and let 


u*j - (Uri) 


Qj> 


'• Q.< 


~ CTk 3 ) , -ft 


fo*- j e j 


(3.79) 


(3.80) 


(3.81) 


(3.82) 


then equation (3.77) will be satisfied. Therefore the 
following equation is used to estimate the value of 


(u kj - c 'k* T kj )/ VV 



if I»j S* 


if !*<• 


(3.83a) 
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where 



\ Tn w] A 


(3.83b) 


An upper limit natural frequency constraint is given by 

J* - U ~ t! i o 

In this case, the optimality criterion, the redesign equa- 
tion, and the estimation of the value [ (U^j - q^*T^ j ) /T^Wj ] 
can be immediately obtained from equations (3.75), (3.76), 

and (3.83), respectively, by replacing the term 

< u kj - W - (u kj - w- 
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CHAPTER IV 


NUMERICAL EXAMPLES 

Computer programs that generate basis design vectors 
were written implementing some of the optimization pro- 
cedures presented in Chapter III. These programs were 
coded in FORTRAN H and were run on an IBM 360/91 compu- 
ter. An optimization program called CONMIN was used 
to obtain the final results reported herein. The 
program CONMIN, developed by Vanderplaats [17J, is 
based on a modified feasible directions method. 

Several design examples are presented here to 
illustrate the effectiveness of the method developed 
in this study. These examples include two and three 
dimensional trusses, and in each example, except for 
the first one, stress, displacement and minimum size 
constraints are included. Some multiple load con- 
dition cases are also considered. 

In order to make the method more effective, an 
approximation technique is employed for estimating 
stresses and displacements during the generation of 
basis design vectors. This technique, which is based 
on using first order Taylor series expansions to expli- 
citly approximate stresses and displacements in terms 
of reciprocal design variables (see Appendix D) , signi- 

PRECEDING PAGE BLANK NOT FILMED 
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ficantly reduces the number of structural analyses needed 
to generate the design basis vectors. Hereafter, the 
optimization method combined with the Taylor series 
approximation technique will be referred to as revised 
method, and the method without the use of the Taylor 
series approximation technique will be designated as 
the ordinary method. The effectiveness of the revised 
method is demonstrated by comparison with the ordinary 
method in some examples. 
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4.1 9 Bar Truss 


The first example problem is a nine bar space truss 
{see Fig. 10) which is studied to demonstrate the appro- 
priateness of the optimization procedure developed in 
Chapter III in comparison with those which were given 
by Gellatly [9] and by Venkaya [11 j (see Appendix E) . 

For the sake of simplicity, only generalized stiffness 
constraints which can be called "total strain energy 
constraints" (see Appendix C) are considered. The 
material properties and the specified value of con- 
straints (upper limit on total strain energy) are given 
in Fig. 10. For this example, two distinct cases are 
considered, and the load conditions for each case are 
given in Table 1(a). Design variable linking is used 
to impose symmetry with respect to both the x-z and 
y-z planes, and the number of design variables is three. 

Results for these two cases are summarized in Table 
1(b). In case 1, the minimum weights obtained are essen- 
tially the same, although Venkaya' s design is heavier 
than the others by 6%. In case 2, however, the design 
obtained by the present method is lighter than the 
others by almost 20%. It is also noted that in the 
present design, both constraints almost reach the spe- 
cified upper limit, but in the other designs the total 
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Material : 

Minimum Size: 

Maximum Size 
Energy Limits: 

Figure 10. 9 Bar Truss 


• 7 

Aluminum, E = 10 psi, p = 0.1 pci 
0.01 in^ 

None 

100 lb-in on both load conditions 
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Table 1 Design Data and Results for Example 4,1 


Oi 


(a) Load Conditions (lb) 



Load 


Direction 

Case 

Condition 

Node 

X 

Y 

Z 


1 

5 

2000.0 

0.0 

-3000.0 

1 


6 

0.0 

O 

• 

o 

-3000.0 


2 

5 

O 

« 

o 

4000.0 

0.0 

— 

f 

j 

6 

0.0 

-4000.0 

O 

• 

o 


1 

1 

5 

3000.0 

0.0 

o 

• 

o 

2 


6 

-3000.0 

0.0 

0.0 


2 

5 

0.0 

4000.0 

0.0 

L_ 


6 

0.0 

-4000.0 

0.0 


(b ) Summary of Results 


Case 

Method 

No. of 
Analyses 

Weight 

(lb) 

Element Size (in^) 

Strain Energy 

1 

2 5 

6 9 

Ld. 1 

Ld. 2 

■ 

Venkaya 

10 

30.5 

0.678 

1.180 

0.158 

100.0 

83.2 

Geliatly 

9 

28.8 

0.736 

0.966 

0.173 

100.0 

92.6 

Present 

10 

28.2 

0.776 

0.883 

0.197 

100.2 

101.7 

2 

Venkaya 

8 

22.8 

0.875 

0.918 

0.01 

30.7 

100.0 

Geliatly 

8 

22.7 

0.869 

0.911 

0.01 

30.8 

100.5 

Present 

8 

19.1 

0.275 

0.911 

0.01 

95.5 

100.5 




strain energy under load condition 1 is only 31% of the 
specified limit, which makes these designs much heavier 
than the present one-. 

As is obvious from Appendix E, both Venkaya’s and 
Gellatly's methods are approximate, and they do not 
guarantee that all the active constraints achieve their 
limits at the same time. It follows that the design 
obtained by these methods cannot be expected to neces- 
sarily be optimal, and in some cases, the results may 
be rather far away from the optimal design . As illus- 
trated by this example, on the other hand, the present 
method is able to overcome this shortcoming and it can 
be expected to produce better results. 
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4.2 10 Bar Truss 


The second example problem is the familiar ten bar 
planar truss (see Fig. 11) for which results have been 
previously reported in [9], [11] and [14]. The material 
properties, stress limits and minimum sizes are given 
in Fig. 11. For this example four distinct cases are 
considered. In case 1-a, the truss is subject to a 
single load condition consisting of 100 Kip downward 
loads applied at joints 2 and 4 (see Fig. 11) and no 
displacement limitations are imposed. Case 1-b is the 
same as case 1-a but with vertical displacement limits 
of ± 2.0 in. imposed at all joints. In case 2-a the 
truss is subject to a single load condition consisting 
of 150 Kip downward loads applied at joints 2 and 4 as 
well as 50 Kip upward loads at joints 1 and 3, and no 
displacement limitations are imposed. Case 2-b is the 
same as case 2-a but with vertical displacement limits 
of ± 2.0 in. imposed at all joints. No design variable 
linking is employed in this example, therefore, the num- 
ber of design variables is ten. 

Results for cases 1-a and 2-a are summarized in 
Table 2, where part (a) of the table contains the results 
obtained by the ordinary method, part (b) of the table 

t 

contains those obtained by the revised method, and part 
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360 


360 



7 

Material : Aluminium, E = 10 psi, p = 0.1 pci 

. 2 

Minimum Size : 0.01 in 

Maximum Size : None 

Load Conditions 

Case 1 : P 1 = 100 kip ' P 2 = 0 

Case 2 : = 150 kip, P2 ~ 50 kip 

Stress Limmits : 25,000 psi 


Fig. 11 10 Bar Truss 
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Table 2 Results for Example 4.2 (1) 


(a) Results obtained by Ordinary Method 


Case 

NO. of 
Analyses 

Weight 

(lb) 

Element Size (in 2 ) 

1 

2 

3 

4 

5 

6 

7 

8 

9 ! 

10 

1-a 

12 

1591.7 

7.88 

0.11 

8.12 

3.89 

0.1 

0.11 

5.82 

5.49 

5.50 

0.16 

2-a 

8 

1664.4 | 

5.94 

0.1 

10.06 

3.95 

0.! 

2.05 

8.56 

2.75 

5.58 

i 

0.1 


(b) Results obtained by Revised Method 


Case 

No. Of 
Analyses 

Weight 

(lb) 

Element Size 

(in 2 ) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1-a 

6 

1593.4 

7.94 

0.1 

8.06 

3.94 

0.1 

0.1 

5.74 

5.57 

5.57 

0.1 

2-a 

6 

1664.6 

; 5.95 

0.1 

10.05 

3.95 

1 0.1 

2.05 

8.56 

i 

2.76 

5.58 

0.1 


(c) Stresses for Results shown in (a) and (b) 



Case 

Stresses (psi) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

(a) 

1-a 

25088. 

19098. 

-24914. 

-25173. 

-871. 

19098. 

24831. 

-25180. 

25173. 

-19098. 

z-a 

25012. 

13097. 

-24993. 

-25006. 

-112. 

24998. 

24989 . 

-25035. 

25006. 

-18522. 

(b) 

1-a 

24991. 

15516. 

-25006. 

-24997. 

76. 

15516. 

25010. 

-24980. 

24997. 

-21943. 

2-a 

24996. 

13094. 

-25002 

-24999. 

38. 

25002. 

25004. 

-24988. 

24999. 

-18454. 






(c) of the table displays the stresses corresponding 
to the designs shown in (a) and <b) . It is apparent 
from these results that the approximate analysis tech- 
nique is rather effective in reducing the total number 
of analyses needed to obtain an optimum design. The 
results shown in part (b) of Table 2 indicate that the 
minimum weights obtained are just the same as those re- 
ported in Ref. [ 14] , and it is interesting to note 
that the fully stressed designs obtained coincide with 
the actual optimal designs in this case. 

Table 3 shows the results of generating basis de- 
sign vectors for displacement constraints. Part (a) of 
the table contains the results obtained by the ordinary 
method, and part (b) of the table contains the results 
obtained by the revised method. Parts (c) and (d) of 
the table display the displacements corresponding to 
the designs shown in (a) and (b) , respectively. The 
results obtained by both methods are essentially the 
same, however the number of analyses required is cut 
in half when the revised method is employed. Observing 
the optimization process, it was found that this ex- 
ample exhibits a particularly interesting behavior. 

This behavior is well illustrated by following the 
iteration history generated by the ordinary method. 

In case 1-b, initially the vertical displacement at 
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Table 3 Results for Example 4.2 (2) 

(a) Basis Design for Displacement Constraints obtained by Ordinary Method 



NO. of 

Weight 




2 

Element Size (in 

) 




Case 

Analyses 

(lb) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1-b 

14 

5044.9 

31.14 

0.1 

22.54 

15.42 

0.1 

0.41 

5.83 

22.07 

21.81 

0.1 

2-b 

14 

4502.3 

24.48 

0.1 

23.49 

14.03 

0.1 

0.97 

9.31 

14.51 

19.85 

0.1 


(b) Basis Design for Displacement Constraints obtained by Revised Method 


No. of 

Weight 


Element Size 

(in ) 

Analyses 

(lb) 

1 

2 3 4 5 

6 7 8 9 10 

8 

5034.4 

31.03 

0.1 22.50 15.32 0.1 

0.79 5.80 21.93 21.67 0.1 


4504.3 24.52 0.1 23.48 14.05 0.1 0.961 9.27(14.54 19 ♦ 87| 0.1 


(c) Displacements for Design (a) 


(d) Displacements for Design {b) 


X 


0.095 


•0.541 


0.233 


■0.309 


0.381 I -0.205 


■1.942 -0.636 


-1.974 


-1.937 


0.092 


■0.550 


0.2 37 


-0.314 


1.996 I -0.391 


- 2.000 - 


0.228 1 -0.617 


0.375 


























joint 2 was assumed to be critical. At the 8th it- 
eration, the vertical displacement at joint 1 be- 
came most critical, and at the 9th iteration, a feas- 
ible design of weight 5227.1 lb. was obtained. However, 
at the 10th iteration, the weight jumped up to 7579.0 
lb. and the vertical displacement at joint 4 also became 
critical. After five more iterations the design weight 
of 5044.9 lb. was obtained, and for this design all 
three of the previously mentioned displacement constraints 
were critical. The iteration history for this case 
(Example 4.2 , case 1-b) is shown in Fig. 12. In 

case 2-b, the vertical displacement at joint 2 was ini- 
tially assumed to be critical. In the 6th iteration, 
the vertical displacement at joint 4 became critical and 
seriously violated. Consequently, the weight suddenly 
increased in the 7th iteration and the displacement at 
joint 2 became most critical again. Finally a design 
weighing 4502.3 lb. was obtained, and for this design 
both the constraints were critical. The iteration history 
for this case (Example 4.2 case 2-b) is shown in 

Fig. 13. A similar phenomenon was reported in Ref. [9] 
for this same problem. The dramatic rise in weight may 
be associated with a major redistribution of the internal 
forces in the structure. 

The final design for cases 1-b and 2-b are summarized 
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Displacement (in) 



Fig. 13 


Iteration History of Example 4.2, Case 2-b 


Weight (lb) 


together with actual stresses and displacements for the 
design in Table 4. In both cases the minimum weights 
obtained are essentially the same as those previously 
reported in Refs. [11] and 114]. The respective criti- 
cal constraints are the tensile stress in element 5 
and the vertical displacement at joints 1 and 2 in case 
1-b, and the tensile stress in element 5 as well as the 
vertical displacement at joint 2 in case 2-b. 
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Table 4 Results for Example 4.2 (3) 


{a) Final Design 


Case 

No. of 
Analyses 

Weight 

(lb) 

2 

Element Size (in ) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1-b 

7 

5077.6 

30.53 

0.13 

23.08 

15.08 

0.13 

0.73 

7.45 

21.56 

21.33 

0.13 

2-b 

7 

4708.6 

23.88 

0.14 

25.32 

13.99 

0.14 

1.96 

12.57 

13.73 

19.78 

0.14 


(b) Stresses for Design (a) (c) Displacements for Design (a) 


Node 

1-b 

2-b 

X 

Y 

X 

Y 

1 

0.195 

-2.000 

-0.031 

-1.099 

2 

-0.546 

-1.992 

-0.609 

-1.999 

3 

0.240 

-0.728 

0.233 

-0.637 

4 

-0.307 

-1.627 

-0.349 , 

-1.528 


Case 

Element 

1-b 

2-b 

1 

6663. 

6469. 

2 

, 

-1257. 

-7334. 

3 

-8519. 

-9698. 

4 

-6642. 

- 7223. 

5 

24962. 

24747. 

6 

-222. 

25000. 

7 

18334. 

16371. 

8 

-6782. 

-5611. 

9 

6642. 

7223. 

10 


10372. 


(psi) 












4.3 2 5 Bar Truss 


The third example problem is a twenty five bar space 
truss (see Fig. 14) for which results have been previ- 
ously reported in Refs. [9], and [ 14 3 • The material proper- 
ties, tensile stress limits, displacement limits and mini- 
mum size constraints are given in Fig. 14. The allowable 
compressive stress limits are listed in Table 5 (b) , and 
they correspond to those given in Ref. 19]* The struc- 
ture is subject to two distinct load conditions as given 
in Table 5(a). Displacement limits of ± 0.35 in. are 
imposed on all joints in all directions. Design variable 
linking is used to impose symmetry with respect to both 
the x-z and y-z planes, and the number of design var- 
iables is eight. 

Table 5 (b) shows the basis design vector for stress 
constraints which was obtained by the ordinary method. 

The basis design vectors for displacement constraints 
are obtained by both the revised and the ordinary methods 
respectively. The results are shown in Table 6 (a) and 
their iteration history is shown in Fig. 15. In both the 
designs shown in Table 6 (a) , critical displacements were 
those in the y direction at joint 1 and 2 under both 
the load conditions. The final design obtained is shown 
in Table 6 (b) , and the weight is very close to the lightest 
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Material 
Stress Limmits 

Displacement Limmits 

Minimum Size 


. . 7 

Aluminium, E=10 psi , p=0.1pci 
Tension 40,000 psi 
(See Table 5 for Compression) 
±0.35 in 

(All Nodes, All Directions) 
0.01 in 2 


Fig. 14 25 Bar Truss 
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Table 5 Results for Example 4.3 (1) 


(a) Load Conditions 


Load 


Direction 

Condition 

Node 

X 

Y 

z 


1 

1000.0 

10000.0 

-5000.0 

1 

2 

0.0 

10000.0 

-5000.0 


3 

500.0 

o 

• 

o 

O 

• 

o 


6 


0.0 

o 

• 

o 

2 

| 

" 

o 

• 

o 

20000.0 

-5000.0 

i 


0.0 

-20000.0 

-5000.0 


(b) Basis Design for Stress Constraints obtained by Ordinary Method 


No. of 
Analyses 

Weight 

(lb) 

2 

Element Size (in ) 

1 

2^5 

6^9 

10, 11 

12, 13 

14^17 

18^21 


4 

344.0 

0.01 

1.250 

1.115 

0.01 

0.01 

0.552 

1.648 

1.336 

Allowable Stress 



-17305. 

-35092. 

-35092. j 


■ - - - ■ 

-6959. ; 

-11082. 


(psi) 



















Table 6 Results for Example 4.3 (2) 


(a) Basis Designs for Displacement Constraints 


Method 

No. of 
Analyses 

Weight 

(lb) 

2 

Element Size (in ) 

1 

2^5 

6^9 

10,11 

12,13 

14^17 

18^21 

22^25 

Rev. 

4 

543.6 

0.01 

2.082 

3.032 

0.01 

0.01 

0.662 

1.656 

2.569 

Ord. 

8 

543.6 

0.01 

2.082 

3.032 

0.01 

0.01 

* 11 

0.662 

1.656 

2.569 


(b) Final Design 


No. of 

Weight 

2 

Element Size (in ) 

Analyses 

(lb) 

1 

2^5 

6^9 

— 

10,11 

12,13 

14^17 

18^21 

22^25 

9 

551.6 

0.010 

2.112 

3.063 

0.010 

0.010 

0.674 

1.690 

2.602 


(c) Final Designs obtained by using Quasi Basis Designs 


Case 

No. of 
Analyses 

Weight 

(lb) 

2 

Element Size (in ) 

1 

2^5 

6^9 

— 

10,11 

12,13 

14VL7 

18^21 

22^25 

1 

‘ 

6 

584.9 

0.013 

2.442 

3.169 

0.013 

0.316 

0.663 

1.947 

2.396 

2 

5 

563.0 




0.013 

0.013 



2.547 
















weight previously reported. The final critical con- 
straints are the stresses in elements 19 and 20 under 
load condition 2, and the displacements in the y direc- 
tion at joints 1 and 2 under load condition 1. 

The designs obtained by using another two sets of 
basis design vectors are shown in Table 6(c). In case 
1 the basis design vectors are obtained by executing 
only one iteration in each procedure for generating the 
basis design vectors previously used. In case 2 the 
basis design vectors were obtained in such a way that 
the redesign procedure was carried out using the Taylor 
series approximations after two exact structural analyses. 
The weights obtained are higher than that of the exact 
design by about 6% and 2%, respectively, however these 
results have practical significance since the computa- 
tional effort required to obtain them is significantly 
less than that needed to produce the results given in 
Table 6 (b) . 
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4.4 72 Bar Truss 


The last example problem is a seventy two bar space 
truss for which results have been previously reported 
in Ref. [9], [11] and [14]. Figure 16 shows the geo- 

metry of the structure, and the node as well as element 
numbering system is illustrated in detail for the upper 
tier. The material properties, stress limits, displace- 
ment limits and minimum sizes are given in the same figure. 
The structure is subject to two distinct load conditions 
as given in Table 7. Displacement limits of ± 0.25 in. 
are imposed on all joints in all directions. Design 
variable linking is employed, and the number of indepen- 
dent design variables is sixteen. 

The basis design vector for stress constraints is 
given in Table 8(a), and that for displacement constraints 
is given in (b) . Both the basis design vectors were ob- 
tained by the ordinary method. The final design obtained 
is listed in Table 9(a). The weight for the design is 
384.8 lb., and is very close to the lightest weight pre- 
viously reported in [11]. The final critical constraints 
are the compressive stresses in elements 1,2,3 and 4 
under load condition 2 and the displacements in both the 
x and y directions at joint 1 under load condition 1. 

The design obtained by using another set of basis design 
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Material 
Stress Linmuts 
Displacement Limmits 

Minimum Size 


Aluminium, E=10psi, p=0.1pci 
±25,000 psi 
±0.25 in 

(All Nodes, All Directions) 
0.1 in 2 


Fig. 16 72 Bar Truss 




Table 7 Load Conditions for Exampl 


Load 

Condition 


Node 


Direction 


X 


5000.0 


0.0 


0.0 


0.0 


5000.0 


0.0 


0.0 


0.0 


0.0 


0.0 











Table 8 Results for Example 4*4 (1) 
(a) Basis Design for Stress Constraints 


Weight 


Element Size (in ) 


I'M 

5^12 

13^16 

h-» 

H* 

00 

19^22 

23^30 

31^34 

35,36 

0.189 

0.1 

0.1 

0.1 

0.190 

0.1 

0.1 

0.1 

37^40 

41M8 

49^52 

53,54 

55^58 

59^66 

67^70 

71,72 

0.199 

0.1 

0.1 

0.1 

0.294 

0.1 

0.1 

0.1 


(b) Basis Design for Displacement Constraints^ 


No, of Weight 


2 

Element Size (in ) 


2 3^-30 

31^34 

35,36 

0.502 

0.1 

0.1 

59-V66 

67^70 

71,72 

0.498 

0.1 

0.1 


























Table 9 Results for Example 4.4 (2) 

(a) Final Design obtained using ordinary method Basis Designs 


f Weight 
Analyses (lb) 


Element Size (in 

•) 


13-16 

17,18 

19-22 

23-30 

31-34 


0.154 0.540 0.410 0.556 0.701 0.510 0.109 0.109 


384.8 


41-4 8 


0.505 


55-58 

1.912 


67-70 71,72 


0.109 


(b) Final Design obtained by using Quasi Basis Design 


1^45 *12 


No. of 

Weight 

Analyses 

(lb) 


35,36 


1.404 0.525 0.347 0.484 1.601 0.488 0.116 0.116 


60.4 


41^48 

0.486 


0.116 2.372 


59^66 67V70 71,72 


0.116 








































vectors is shown in Table 9 (b) . These basis design vec- 
tors are obtained by executing only one iteration in each 
procedure for generating the basis design vectors previ- 
ously used. The weight obtianed is 460 lb. and it is 
heavier than the other design by almost 20%. 

Finally the number of structural analyses and CPU 
runtime are summarized in Table 10 for each example 
problem. 



Table 10 Summary of Number of Analyses and Run Times 
























CHAPTER V 


CONCLUSIONS 

In this study/ the primary effort was focused on 
the derivation of optimality criteria and the develop- 
ment of rational recursive redesign procedures based 
on the optimality criteria. Optimality criteria, which 
are equivalent to the necessary conditions for local 
optimality, were derived for a general case, so that 
they could easily be specialized for each type of 
behavior constraint considered. Special attention 
has been given to developing design optimization procedures 
for basis vector generation that: (1) are rational; (2) 

are efficient; (3) yield feasible basis design vectors. 

The recursive redesign procedures presented are based 
on optimality criteria concepts. Furthermore the pro- 
cedures presented are rational and they do not exhibit 
the shortcoming present in some of the previously reported 
methods (for example, see [9] and [llj), for cases in- 
volving multiple constraints (of the same behavior type) 
and multiple load conditions . As shown in example problem 
4.1, this was accomplished by the use of a minimum square 
method to estimate the optimal values of the multipliers. 

It is also noted that the estimation procedure requires 
little additional effort because the values of the mul- 

RECEDING PAGE BLA NK NOT FILMED 
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tipliers can be obtained by an equivalent linear program- 
ming method. The procedures presented for generating 
basis design vectors are efficient. To attain this 
efficiency, a method for deleting noncritical or less 
critical constraints was adopted. For the example prob- 
lems, the constraint deletion technique worked well be- 
cause the number of active constraints at the optimum 
design is quite small in comparison with the total 
number of constraints, and most of the active constraints 
remain active during the entire design process. Further- 
more, first order Taylor series expansions with respect 
to reciprocal design variables were used to provide 
explicit approximate representations for stresses and 
displacements. This high quality approximation was 
very effective in reducing the number of actual analyses 
needed to obtain an optimum design. Finally it should 
be noted that care has been taken to avoid the genera- 
tion of infeasible basis design vectors such as the 
one shown in Fig. 7. This was achieved by imposing 
additional limitations on the values of the multipliers 
corresponding to violated constraints. 

With regard to the results for the example problems, 
one notable thing is the accuracy of the approximate 
optimum designs. The obtained weights corresponding 
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to the final designs obtained are very close to (less 
than 1% above) the best previously reported results for 
these example problems. It is also noted that these 
designs were obtained using only two basis design vec- 
tors in each case. The results reported here suggest 
that basis design vectors generated by optimality 
criteria and stress ratio methods may frequently span 
a subspace such that the application of mathematical 
programming methods to the reduced problem, cast in 
terms of a few generalized design variables, provides 
an efficient hybrid method for obtaining an excellent 
upper bound approximation of the optimum design. The 
trade off between the effort expended to refine the basis 
vectors (i.e., converge the subproblems) and their qual- 
ity, with respect to spanning a subspace containing a 
good upper bound approximation of the optimum design, 
is an open question which will require further study. 

This trade off must be investigated thoroughly before 
hybrid methods of structural optimization can become 
efficient tools for the design of large practical struc- 
tural systems subject to a wide range of behavioral con- 
straint types. It should be emphasized that the signifi- 
cance of the results reported here is that they establish 
the feasibility of the hybrid method concept and they 
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indicate the promise that these methods hold for achieving 
high efficiency. 

Based on the study presented here, the following 
conclusions have been reached: 

(1) The hybrid method appears to be one of the most 
efficient methods especially for large scale struc- 
tural problems. 

(2) Basis design vectors for hybrid methods can be gen- 
erated efficiently by the optimality criteria methods 

(3) The stress ratio method appears to be adequate for 
generating the design basis vector containing infor- 
mation relative to the stress constraints. 

(4) In many cases the optimum designs obtained for each 
subproblem form a reduced basis that spans a sub- 
space containing a good upper bound approximation 
of the optimum design. 

(5) Taylor series expansion with respect to linked reci- 
procal design variables can be used in the context 
of optimality criteria and stress ratio methods and 
they produce considerable improvement in efficiency 
by reducing the number of actual analyses needed 

to achieve convergence. 

As a result of this study, it is suggested that the 
following additional work may be of interest: 



( 1 ) 


The method developed for buckling and natural fre- 
quency constraints should be implemented and nu- 
merical results should then be generated and examined 

(2) A method to prevent convergence to nonoptimal points 
such as the 'one shown in Fig. 8 should be sought. 

(3) A trade off study employing various strategies that 
combine optimality criteria methods, or stress 
ratio methods, with the Taylor series expansion 
technique for approximate analysis should be car- 
ried out. The trade-off will be between the effort 
expended to obtain the basis vectors and their qual- 
ity, with respect to spanning a subspace containing 
a good approximation of the optimum design. 

and finally 

(4) The method should be extended to other types of 
behavior constraints as well as a broader class of 
structures. 
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APPENDIX A 


FARKAS ’ LEMMA AND DERIVATION 
OF EQUATION (3.17) 


Farkas 1 Leiruna Ref. 18 is given as follows: 

Let {? Q , k e an arbitrary set of vectors. 

There exists 8^ ^ 0 such that 

- S, -> , 

Po “ Z— Pi 

if and only if 



for all y satisfying 

__ 

Pi z o , i » i. 2, r 


(A.l) 


From the lemma/ the following relation can easily 
be derived by replacing P by (-P^) , i = l,2,...,r, 
if and only if 


P. T 7 S 0 

for all y satisfying 

P^ 7 1 0 / 1 * 1.2,--*/ Y" 


there exists 8^ >_ 0 such that 


2L/ 6; Pi = « 

A*' ' 


t (A. 2) 


) 
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A set of equations (3.16) can be rewritten in the following 
matrix form 

SW =■ sd 

for all SD satisfying 

sf* = Vj>- SD 4 0 , for I? 6 K«t 
- SD- » I* - $D id , fov 3 ■€ vJmSM 
S D j ** j \T »«yt 

where 

AfJ ~ ( 'Uj ~ i , ■1a /' ■ • * * , 'US' i mov ) 

s'? - ( SDi, SD 2/ ■ ■ , SDudv > 

= ^ Jm * J fc*/ *** - J^npv ^ 

and Jj are unit vectors in which only the jth element 

has the value of -1 and 1, respectively, and all the other 
elements are zero. 

Using the relation (A. 2) and replacing P q by w, 
y by SD, P i by Vg k , l j and J , 3 ± by ^ y. and ru , 
it is proved that the following relation is mathematically 
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equivalent to that given by (A. 3), 


* 2 _ *]* + . 1 - ft % - 2 _ 1 ; 


= 0 


(A. 4) 


where y. and y . are nonnegative. 

*3 3 

In order to include inactive constraints into the relation 
(A. 4), we introduce the null multipliers such that 


> jt = 0 , -for fc Kact 


Using the null multipliers, the relation (A. 4) can be 
rewritten as 


mc 


W + TL Xk *3* 

fr-» 0 


\2-*i ^ 

'['to mi* 


■ ^ *li Ji 


= o 


(A. 5) 


From equation (A. 5), we get 
Me 


+ I- J*,; = 0 ' f or 1 
0 J 


£ J 


NC 




i + }r x K = /V 1 * J " 


H C 


+ Jki *-*ij ' j ^ J| 






4 


*KA# 


Noting that yu and are nonnegative, we get the follow- 
ing relations 
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= o , f*r j « T 


A <i + ^ a ~ 0 ' ^ i * J * 


f oy ^ ^ j , 


where 


* 0 , for fr * K a ct 

* O , for -<C K att 


This relation is just the same as that given by (3.17) 
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APPENDIX B 


MATHEMATICAL CONSIDERATION 
ON THE HYBRID METHOD 

Here some mathematical considerations for the hybrid 
method (basis design vector method) are presented. As 
shown in Chapter II, the structural optimization problem 
considered in this paper is given in the following form 


Minimize W = W(D) 

subject to (B.l) 

4 0 , fe - 1. 2, K 

where D represents an M dimensional design variable vector, 
and g^(D) includes both behavior and size constraints and 
K denotes the total number of constraints. 

Let D^, D 2 ,...,D n , N < M, be an arbitrary set of 
M dimensional vectors, and define a set of new design 
variables 0 ^ such that 




(B . 2 ) 


Assume that we get the optimal solution of problem (B.l) 
for % , let it be then the following optimality criteria 

must be satisfied 

e 0 , 1 * l. 2. / M (B . 3) 

^ 0 ; i. i F -i A* 1 


tr« I 


5 © 


where 


N 


D 


fr «i 


(B.4) 


103 



From equation (B.2) 


-i_ = Y TV* — ' - ■ 

Y *V 


(B. 5) 


where denotes the ith element of D^. From equations 

(B.3) and (B.5) 




(B.6) 


If D* is the actual optimal solution of problem (B.l) , 
the term in the parenthesis should be zero for all i. 
However, equation (B.5) does not guarantee it because 
the number of equations N is less than the number of the 
terms in the parentheses which is M. It follows that 
the actual optimal solution does not necessarily exist 
in the subspace defined by equation (B.2) for any choice 
of Dj , j = 1,2,...,N. Therefore, the solution given by 
equation (B.4) must be an approximation of the actual 
optimal solution. 

Next we will consider the case that Dj is given 
by the optimal solution of the jth subproblem which is 
defined as follows 
Minimize W - W(D) 
subject to 

i 0 , k - 1,2,-, Kj 

C P ) 4 0 / k i , 2, -* - , K o 
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where g^fD) represents the constraints which are to be 
imposed only in the jth subproblem, and g o ^(D) denotes 
the constraints to be imposed in every subproblem. Let 
Kj and K q be the number of corresponding constraints, 
respectively, and investigate the conditions under which 
the method will give the actual optimal solution. From 
the assumption, the following optimality criteria must 


be satisfied for each j , 


LimSLI 

7>I>i 


Wk 

k'1 1 * 


i = 1,2, . . . ,M (B . 7) 

If D* defined by equation (B.4) is the actual optimal 
solution, then 

+ TL hr ♦ Z fr =0, 

A « I, 2, M, K“ ilK.' •+ 

1 

If equation (B.8) can be expressed by the linear combina- 
tion of equation (B.7), there exist A^ and satisfying 
equation (B.8), and consequently 6* gives the actual 


optimal solution. Generally, this is possible only 
when g^ is a linear function of D for all k, and W is 

either of the following 

v 

(1) linear function of D 

(2) W = Is6 T [A]B + W q 
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where [A] is the MxM matrix and W q is constant. Under 
these conditions the hybrid method may give the actual 
optimal solution. In other cases, it may depend on the 
problem itself, and general conditions have not yet been 
derived. 
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APPENDIX C 


OPTIMALITY CRITERIA FOR 
GENERALIZED STIFFNESS CONSTRAINTS 
Generally speaking, the total strain energy stored 
in a structure represents an inverse measure of the 
stiffness of the structure. Therefore, the stiffness 
requirement can be set up by restricting the value 
of total strain energy. This is called "generalized 
stiffness constraint" [11], and it is usually given 
by 

<\ - \ Ct T CK3 u - U* i 0 (c.i) 

where U* denotes the specified upper limit of total 
strain energy. 

Differentiating both sides of equation (C.I) with 
respect to A^, we get 

ffc - (C.2, 

Substituting equations (3.39) and (3.43) into equation (C.2) 
yields 

ffc - -il’Wt 

Using the relation given by (3.36), we get 

= - 4 - 


(C. 3) 


(C.4) 
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From equation (3.17), we obtain the optimality criteria 
for generalized stiffness constraints such that 


" "T £ 

** I 


S 


u 


T [)?i3 u « o, -for j -e j (c.5) 


Multiplying both sides of equations (C.5) by D ^ , we 
get the following standard form, 

W; - > U j =0, -f or j < O’ ( c - 6 > 

whe re 

Uj - - Z. w T De;] u 

(I 

then Uj denotes the total strain energy stored in the 
jth group of elements. 

For multiple constraints, the criteria can be gen- 
eralized as 


Me 


W; - X* Ufej = 0, -for j € J 


(C . 7) 




where 




2 0/ ■f' >y fe ”€r act 

~ 0 , -f 6y ^ ^ Katt 


and U. . represents U. for the kth constraint. 
K J J 
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The redesign equation can be obtained in 
similar to that used for displacement constraints, and it is 


< t>i >*» = C .i > 


d 


'd 


where 


l (C . 8) 


ML 


Cj -ItH 


U 


1 . 
2 




W, 


< 

If size constraints are imposed, equation (3.20) must be 

used together with equation (C.8). 

The value of U. ./W. can be estimated in the following 
k j 3 

manner. For active constraints 




uf 


mpv 

i" 


(C . 9a) 


(C. 9b) 


where (U ) represents the total strain energy of a 
k s 

whole structure at the sth iteration. 

Let 


u 


u* 


*j “ CUk), 


(C.10) 


then equation (C.9a) is satisfied. Thus we will use the 
following equation to estimate the value of U^/W. 
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s 


(C.ll) 




Uk 

CU*>, 
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APPENDIX D 


LINEAR APPROXIMATION OF 
STRESS AND DISPLACEMENT 

It has been recognized that structural behaviors 
such as stress, displacement and so on can be estimated 
by using a first order Taylor series expansion, which 
is given in the following form 

■je I « A | V 

where f(x) is an arbitrary differentiable function of 
variable x = {x^}, and x q = {Xj Q } an ar bitrary given 
point. Applying equation (D.l) to displacement yields 


N 

uO?.) + 21 


am 50 
axj 


C X j Xj o ) 


(D.2) 


where x represents an appropriate design variable vector. 
Equation (d. 2) must be applicable to any displacement, 
thus we get 



uOO + 

i-' 





(D.3) 


It has also been found that the use of reciprocals 
of the sizing type design variables is very effective 
in increasing the accuracy of this estimation. There- 
fore, Xj is selected as 
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*i - -57 ' i - i-'*— » 

From equation (D.4) 

^ 1A ^ ^ tt 

>Xj <* 


(D.4) 


(D.5) 


Substituting equations (D.4) and (D.5) into equation (D.3) 
yields 

_ y -?* SUCgp) f J L\ (D.6) 

P T >E>i V D i T>i- J 

From equations (3.36), (3.39) and (3.43) 


f^r = - [K]" Z_, a.j [tec] w 


» 


(D. 7) 


'i 

Substituting equation (D.7) into (D.6), we get 

u(5) UCD„) 

+ LKiSj ( Z. H S - jr ) (D. 8) 

N j«» <*j ' J 4 7 

If we know u and [K] for design D q , then displacement for 
a new design D can be estimated by using equation (D.8). 

In the displacement method of analysis, stress is 
readily expressed as a function of displacements. There- 
fore stress can also be estimated using equation (D.8). 
For an axial force element, stress is given by 

<f - [S 3 U (D. 9) 


where [S] is a geometrically determined matrix. Sub- 
stituting equation (D.8) into equation (D.9), we get 
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? ( L> ) = 7 c ) 

i- M [KC L S' ^(dj " i ~) 


(D. 10) 
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APPENDIX E 


REDESIGN PROCEDURES BASED ON 
OPTIMALITY CRITERIA PREVIOUSLY PRESENTED 

Here a brief explanation of two representative re- 
design procedures based on discretized optimality criteria 
is presented. Among these two procedures, one was given 
by Venkaya .[11], and the other was given by Gellatly [9]. 

Venkaya has presented a redesign equation for gen- 
eralized stiffness constraints under multiple load con- 
ditions in Ref. [11], which is 

f ± ..u>/ V 

(E ' 1) 

where ou is the ith relative design variable, and A 

(k) 1 

is a scaling parameter. u i is the tota l strain energy 

of the ith element under the kth load condition for the 
relative design, and 1 denotes the weight of the ith 
element for the relative design. s represents the cycle 
of iteration. The weighting parameter c^ is given by 

tfr = ' k * I' 2 ’"*' ? (E . 2) 

where W is the current total weight of the structure and 
is the specified value for the kth constraint, p is 
the number of constraints. 

Equation (E.l) can be rewritten as 


fRECED »G PAgb 


Br, 


■ank 


Not 
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(E.3) 



( k ) 

where is the strain energy of the ith element 

under the kth load condition, and is the weight of 
the element. Under design variable linking, equation 
(E.3) can be modified by using the notation defined in 
this paper as 

w/ ) s (D.4 

J L- 

where U^j represents the total strain energy of the jth 
'group elements , and is the total weight of the group. 

No redesign equation for generalized stiffness con- 
straints has been given by Gellatly. However, he has 
presented a redesign procedure for a combination of 
stress and displacement constraints. The basic concept 
of his method can be summarized as follows: Compute a 

new value of each design variable for each constraint and 
select the largest one for each design variable [9]. 

This concept was applied to the first example problem 
in Chapter IV as Gellatly' s redesign procedure. 
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